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We report on an undergraduate student project initiated in the summer semester of 2004 with the 
aim to establish equations of state for white dwarfs and neutron stars for computing mass-radius 
relations as well as corresponding maximum masses. First, white dwarfs are described by a Fermi 
gas model of degenerate electrons and neutrons and effects from general relativity are examined. 
For neutron star matter, the influence of a finite fraction of protons and electrons and of strong 
nucleon-nucleon interactions are studied. The nucleon-nucleon interactions are introduced within 
a Hartree-Fock scheme using a Skyrme-type interaction. Finally, masses and radii of neutron stars 
are computed for given central pressure. 
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I. INTRODUCTION 



Compact stars, i.e. white dwarfs and neutron stars, are the final stages in the evolution of ordinary stars. After 
hydrogen and helium burning, the helium reservoir in the star's core is burned up to carbon and oxygen. The nuclear 
processes in the core will stop and the star's temperature will decrease. Consequently, the star shrinks and the 
pressure in the core increases. 

As long as the star's mass is above a certain value, it will be able to initialise new fusion processes to heavier 
elements. The smaller the mass, the more the core has to be contracted to produce the required heat for the next 
burning process to start. If the star's initial mass is below eight solar masses, the gravitational pressure is too weak to 
reach the required density and temperature to initialise carbon fusion. Due to the high core temperature, the outer 
regions of the star swell and are blown away by stellar winds. What remains is a compact core mainly composed of 
carbon and oxygen. The interior consists of a degenerate electron gas which, as we will discuss, is responsible for the 
intrinsic high pressure. The compact remnant starts to emit his thermal energy and a white dwarf is born. 

As there are no fusion processes taking place in the interior of the white dwarf, a different force than thermal 
pressure is needed to keep the star in hydrostatic equilibrium. Due to the Pauli principle, two fermions cannot occupy 
the same quantum state, hence, more than two electrons (with different spin) cannot occupy the same phase space. 
With increasing density, the electrons fill up the phase space from the lowest energy state, i.e. the smallest momentum. 
Consequently, the remaining electrons sit in physical states with increasing momentum. The resulting large velocity 
leads to an adequate pressure of the electrons — the degeneracy pressure — counterbalancing gravity's pull and 
stabilising the white dwarf. 

The electrons are the first particles to become degenerate in dense matter due to their small mass. White dwarfs 
consist also of carbon and oxygen, but the contribution of nuclei to the pressure is negligible for the density region of 
interest here. The gravitational force increases with mass. Consequently, a larger degeneracy pressure, i.e. a greater 
density, is needed for stability for more massive white dwarfs. There is a mass limit for white dwarfs beyond which 
even the degenerate electron gas cannot prevent the star from collapsing. This mass limit is around 1.4 solar masses 
— the famous Chandrasekhar mass. For greater masses, the white dwarf collapses to a neutron star or a black hole. 
The radii of white dwarfs are in the range of 10,000 km — about the size of our planet earth. 

Similarly, a neutron star is stabilised by the degeneracy pressure also. The difference is that for neutron stars the 
degeneracy pressure originates mainly from neutrons, not from electrons. 

If the initial mass of an ordinary star exceeds eight solar masses, carbon and oxygen burning starts in the core. 
Around the core is a layer of helium and a layer of hydrogen, both taking part in fusion processes. As long as 
they produce the required temperature to keep the star in hydrostatic equilibrium, new burning processes to heavier 
elements will take place in the core while the fusion of lighter elements will continue in the outer shells. Once iron 
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is produced in the core, there is no burning process left to generate energy as the fusion to heavier elements requires 
energy to be put in. So the mass of the iron core increase while its radius decrease as there is less thermal pressure to 
work against gravity. At some critical mass, the iron core collapses. The density increases so much that the electrons 
are captured by nuclei and combine with protons to form neutrons. Therefore, the atomic nuclei become more and 
more neutron-rich with increasing density. At some critical density, the nuclei are not able to bind all the neutrons 
any more, neutrons start to drop out and form a neutron liquid which surrounds the atomic nuclei. For even larger 
density there will be just a dense and incompressible core of neutrons with a small fraction of electrons and protons. 
Neutrons in neutron stars are then degenerate like electrons in white dwarfs. 

There will be an outgoing shock wave generated by the proto-neutron star due to the high incompressibility of 
neutron star matter. The falling outer layers of the ordinary star bounce and move outwards interacting with the other 
still collapsing layers and generate an overall outward expansion — a core-collapse supernova arises. The neutron 
rich remnant of the core collapse will become a neutron star as long as his mass is less than about two to three solar 
masses, otherwise it proceeds to shrink and becomes finally a black hole. The radii of neutron stars are typically 10 
km — the size of the city of Frankfurt. 

During the summer semester of 2004, the two last coauthors organised a student project about compact stars on 
which we will give an extended report in the following. We derive the equations of state (EoS) for white dwarfs and 
neutron stars, respectively, where we will cover the Newtonian and the relativistic cases for both kinds of compact 
stars. These EoSs are used to determine numerically the masses and radii for a given central pressure. In the first 
chapter we give a short description of the properties of white dwarfs and neutron stars. We then proceed in the 
second chapter with calculations of the EoS for white dwarfs. We assume a Fermi gas model of degenerate electrons 
and outline two different ways which we employed for deriving the equations of state, repeat the basic structure 
equations for stars as well as the corrections from general relativity (the so called Tolman-Oppenheimer-Volkoff 
(TOV) equation) and discuss the Chandrasekhar mass, which is the famous mass limit for white dwarfs. From the 
third chapter on we will concentrate on neutron stars. Again we will have two different approaches for calculating 
the equation of state for neutron stars. After describing a pure neutron star, we will discuss the effect of the presence 
of protons and electrons which will appear in /3-equilibrium. Furthermore, we will introduce a simple effective model 
for incorporating the nucleon-nucleon interactions. We will compare our results to the recent work by Reddy and 
Silbar .lj who are using different numerical methods and another model for the nucleon-nucleon interactions. We 
investigate whether the two model descriptions satisfy causality. Finally, we summarise our findings and give a short 
outline of modern developments in neutron star physics. 



II. WHITE DWARFS 



A. Definitions 



1. Structure equations 



There are two forces acting on the star, one of them is gravitation and the second one arises from the pressure 
(thermal pressure for normal stars and degeneracy pressure for white dwarfs and neutron stars). The pressure and 
force are related in the standard way: 

dF dF 

dp = T = 4^ ' (1) 

where p denotes the pressure, F the force, and A is the area where the force takes effect. For the gravitational force 
one has for a spherical symmetric system: 

dF = - Gdm r2 m(r) , (2) 

with r being the radial distance of a spherical star and 

dm = p(r)dV = p(r)47rr 2 dr . (3) 
Equating both expressions gives the first basic structure equation for stars in general: 

dp = Gp(r)m(r) 

dr r 2 { ' 
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Here G is Newton's Gravitational constant, p is the mass density, m is the mass up to radius r and V is the volume. 
Again using eq. 10, one has also the equation: 

dm , . , 

— — = ptrUirr^ . (5) 
dr 

We can rewrite the mass density in terms of the energy density e: 

e(r) 

P(r) = , (6) 

where c is the speed of light. With these definitions one obtains structure equations which describe how the mass and 
the pressure of a star change with radius: 

dm 4irr 2 e(r) 

dr c 2 

dp Ge(r)m(r) 

T = 2~2 ■ 

dr c^r* 

Here one has to solve two coupled differential equations. Note that while dm/dr is positive, dp/dr must always be 
negative. Starting with certain positive values for m and p in a small central region of the star the mass will increase 
while the pressure decreases eventually reaching zero. One sets m(r = 0) = and, in addition, one has to specify 
some "initial" central pressure p(r = 0) = po in order to solve for eqs. and JSJ. The behaviour of m and p as a 
function of the radius will become important for our numerical computations. 

2. Equations of state 

From the last section, we arrived at two coupled differential equations for the mass and the pressure in a star 
where both depend on the energy density e. One needs to fix now the relation p(e) between the pressure and the 
energy density in addition. Matter in white dwarfs can be treated as an ideal Fermi gas of degenerate electrons. For 
electrically neutral matter an adequate amount of protons to the electron gas has to be added. Furthermore, the 
protons combine with neutrons and form nuclei which will be also taken into account. 

The distribution of fermions as an ideal gas in equilibrium in dependence of their energy is described by the 
Fermi-Dirac-Statistic: 

exp[(E - fij/k B T\ + 1 

where E is the energy, /i is the chemical potential, ks is the Boltzmann constant and T is the temperature. In kinetic 
theory we find the following correlation between the distribution function / and the number density rij = dNi/dV in 
phase space for a certain particle species i: 

= 9 f no) 

d 3 k (27Tft) 3 1 ' 1 ' 

(2iT'h) 3 ' is the "unit" volume of a cell in phase space and g is the number of states of a particle with a given value of 
momentum k. For electrons g equals 2. The number density of species i is given by: 



[2rTt) 

For degenerate fermions the temperature can be set to zero (as (/i — m e c 2 ) 3> fcsT), accordingly ii/ksT goes to 
infinity. To be added to the system, a particle must have an energy equal to the Fermi energy of the system as all the 
other energy levels are already filled. This description matches exactly the definition of the chemical potential, and 
so we can set here /i= Ep. With all these assumptions we can write the distribution function as: 

m = \l t lt E z\ ■ (12) 



for E > E F 

If we ignore all electrostatic interactions, we can write for the number density of the degenerate electron gas: 

- dAk = / k ' 2dk = — rr ■ ( 13 ) 



(27rft)3 (27rft)3 Jo 3ir 2 h 
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As the star should be electrically neutral, each electron is neutralised by a proton. If we assume that the white dwarf 
is predominantly made of 12 C or 16 O, A/Z = 2. The proton and neutrons masses are both much larger than the 
electron's one, so when we calculate the mass density in the white dwarf we can neglect the mass of the electrons and 
just concentrate on the nucleon mass mjv. The mass density in terms of m^r is: 



A 

p = n- m N ■ — , 

where n = n e . With this equation and the dependence of number density on k F we find: 

'3ir 2 pZ" 



k, 



niN A 



(14) 



(15) 



As the momentum of the nucleons is negligible compared to their rest mass, the nucleons do not give a significant 
contribution to the pressure at zero temperature. On the other hand, the electrons behave as a degenerate gas, and 
have large velocities which give the dominant contribution to the pressure. The energy density can be divided in two 
components, one term coming from nucleons and the other one coming from electrons, whereas the contribution from 
nucleons is dominating. The complete energy density can be written as 0: 



e = nm N ^c 2 + e clcc (k F ) . 
where e e \ oc (kF) is the energy density of electrons. With the electron energy 



we can write e e ; ec in the following way: 

£elec(/CF) 



E(k) = yjk 2 c 2 + m 2 c 4 



E(k)k 2 dk 



(16) 



(17) 



8vr 


(2tt^) 3 


8 


n 


(27Tft) 3 











eo 
8 


(2x 



(P^ + mlcy^Pdk 
\u 2 + l)VVdu 



with 



and 



mfc 5 



TT 2 h 3 

x = kp/m e c. 



(18) 

(19) 
(20) 



The factor eo carries the dimension of an energy density (dyne/cm 2 ). The pressure of a system with an isotropic 
distribution of momenta is given by: 



1 8tt 



3 {2-Khf 



kvk dk 



(21) 



where the velocity v = kc 2 /E and the factor 1/3 comes from isotropy. For the electrons we get: 



p(k F ) 



1 


8tt 


r kp k 2 c 2 


3 (2tt?i) 3 


lo E{k) 


8tt 


/ (k 2 c 2 - 




3(2^7i) 3 j 


eo 
3 J 


rkp /m e c 

1 (u 2 + l) 




eo 
24 


"(2x 3 


- 3as)(l + 



k 2 dk 



(22) 
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The energy density is dominated by the mass density of nucleons while the electrons contribute to most of the pressure. 

We want to arrive at an equation of the form p = p(e). In the following, let us consider the two extreme cases: 
i<Cl and x ^> 1, i.e. kp <C m e c and kp 3> w e c, respectively. In the former case, the kinetic energy is much smaller 
than the rest mass of the electrons, i.e. that is the non-relativistic case. Considering the equation for the pressure for 
kp <C m e c one finds that 

PikF) = jJ +1) I udu^-J^ u d u = -{—) = (23) 

With the assumption e = p ■ c 2 one arrives at the following equation of state (EoS) in the non-relativistic limit: 

p « ^ non _ rcl e 5 / 3 (24) 



_ K ( 2,tt 2 Z 



with 

5/3 

157T 2 TO e \mArC 2 A / 

For the relativistic case (&f 3> m e ) we arrive at: 

j»(e) « if rc ie 4/3 , (26) 

where 

y , (2r) 

The relation of the form 

p = ife 7 (28) 

is called a "polytropic" equation of state. With the two EoSs for the relativistic and for the non-relativistic case 
we are in the position to calculate the mass and the radius of a white dwarf with a given central energy density by 
inserting these relations into the structure equations J7J and JSJ. A discussion on general forms of polytropic EoSs 
can be found in |3|]. 

Besides the non-relativistic and the relativistic polytropic EoS there exists also a third case, which becomes impor- 
tant for later considerations the ultra relativistic EoS. Here one assumes for simplicity a relativistic pure electron 
Fermi gas with the total energy density: 

e = e elec (k F ). (29) 
Just like in eq. (fTHf> one can calculate the energy density using 

e(fejr) = e / {u 2 + l) 1/2 u 2 du, (30) 



getting for the relativistic case 



From eq. (|2*2*|l 



one finds for kp 3> m e : 



e = ^ . (31) 

P(k F ) = 4 / (" 2 + 1)~ 1/2 ^ (32) 

P=^ 4 - ( 33 ) 
Therefore, one arrives at the general solution for an ultra-relativistic gas: 

V=\t ■ (34) 
This relation is also very important for neutron stars as a limiting case of their relativistic EoS. 
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3. Chandrasekhar mass 



We will derive in this section the maximum mass for white dwarfs for a polytropic EoS as found by Chandrasekhar 
in 1931 PJ, following Note, that already in 1932 Landau found with simple arguments a mass limit for white 

dwarfs and neutron stars [f| (see also the discussion in e.g. @])- 

By straight-forward algebra the structure equations J7J) and © can be combined to give: 

1 d (r 2 dp\ _^ Gp 



r- 



2 dr V p dr 



The parameter 7 in the polytropic equation is usually rewritten as 7 = 1 + — where n is the polytropic index. 
Given the equation of state p = p(p) in the polytropic form (see eq. (J2HJ) 

p = Ke> = Kp 1 ^ (36) 

one can obtain p(r) by solving eq. (|35|l with the initial conditions: 

p(r = 0) = po jt (37) 

and 

dp 



dr 



= (38) 

(r=0) 



which follows from the condition m{r = 0) = and eq. |JS}. Equation l-iol) can be transformed into a dimensionless 
form with the following substitutions: 

// , -i\jv- (l-n)/n\ 1 / 2 

r = ai with ffl = ^ c(" +1 )/" (39) 

and 

p{r) = p Q n with = 6{r) . (40) 

Here 9 and £ correspond to a dimensionless density and radius, respectively, and a is a constant scale factor. With 
these definitions eq. (|35|l translates to: 

Equation l|41|) is called the Lane-Emden equation for a given polytropic index n. From eq. (|40|l one realizes that 

6(r = 0) = 1. (42) 

Furthermore, with eq. I|38ll the density distribution at the center of the star obeys: 

d6 



dr 



= 0. (43) 

(r=0) 



With the boundary conditions, eqs. I|43|l and l|42() . one can integrate eq. (|41() numerically. One finds for n < 5 and 
7 > 6/5, respectively, that the solutions decrease monotonically with radius and have a zero for a £ = £1, i.e. 0(£i) = 
or p(ri — a£i) = 0, respectively. Hence, the radius of the star is given by R = r\ = a£i which translates to 

R . ( ( ° + tf" , ")' ; '^'" ,fr (44) 
Using eqs. (|4UH and (|39|) the total mass M of the star can be calculated as follows: 

f R 

M = / Anr 2 pdr 



= 47ra> / ¥0 n d£ 

3/2 



= 4^«+3)/n (^^) Pr n)/2n e?|^(6)|. (45) 
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For po we can use eq. (I44|l and get: 



M = And 



(2n+2)/(n-l) /" ( n + jOj^ 



4ttG 



n/(n-l) 



(n-3)/(l-n) t 2|/j/ 



-„)/(!-„) 



The interesting case is the high density limit, i.e. the relativistic case with 7 = |. For this case it follows 

4 
3 
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n=3, 6= 6.89685, | = 2.01824 



(46) 



(47) 



For K we take eq. <|27[) and set 77 = AjZ. Using these values in eq. (|46|l . we get the following expressions for the mass 
and the radius in the relativistic case (see @). For the radius: 



1 , / h 3/2 

R= -(3^)^(6.89685) 



Pcrit 

, Po 



1/3 



with 



and for the mass: 



1 



Pcrit 



M = ^(3vr) 1/2 (2. 01824) ( 



rriNr]m e 3 c 



he 



G 



3/2 



1 



m N r] 



which becomes after plugging in all the values Q, 



M = 1.4312 



(48) 



(49) 



(50) 



(51) 



Here and for the later calculations we chose for tun the neutron mass. For a white dwarf which consists mainly of 12 C 
it is better to take the atomic mass unit (again from Q), which gives a maximum mass of 1.4559 solar masses. It is 
remarkable that in the relativistic case the mass does not depend on the radius R or the central density po anymore. 
As the density in the white dwarf increases, the electrons become more relativistic and the mass approaches the value 
of eq. (|5 II) which is the famous Chandrasekhar mass M c f L . It represents the maximum possible mass for a white dwarf 
for a purely relativistic polytropic EoS. 



4- General Relativity Corrections 



If the stars are very compact one has to take into account effects from general relativity, like the curvature of 
space-time. These effects become important when the factor 2GM / c 2 R approaches unity. We then have to describe 
a compact star by using Einstein's equation: 

G ijlv = ~^~T^ V . (52) 

For an isotropic, general relativistic, static, ideal fluid sphere in hydrostatic equilibrium one arrives at the Tolman- 
Oppenheimer-Volkoff (TOV) equation (see e.g. 0, IE IE m, f° r detailed derivation): 



dp 
dr 



Ge(r)m(r) 



1 



P(r) 
e(r) 



1 



4irr 3 p(r) 
m(r)c 2 



1 - 



2Gm(r) 



1 -1 



(53) 



This equation is similar to the differential equation for the pressure, eq. ©, but has three correction factors. All 
three correction factors are larger than one, i.e. they strengthen the term from Newtonian gravity. The TOV equation 
contains also the factor 2GM/c 2 R which determines whether one has to take into account general relativity or not. 
The corresponding critical radius 

2GM r*^ 
R = — 5— (54) 

<r 

is the so called Schwarzschild radius. For a star with a mass of 1 Mq, we get R » 3km. Hence effects from general 
relativity will become quite important for neutron stars as they have radii around 10 km. The corrections are smaller 
for white dwarfs with radii of R(M = 1M Q ) w 10 4 — 10 3 km but will be checked numerically later. For completeness, 
note, that an anisotropy in the pressure can generate an additional term to the TOV equations even for spherical 
symmetry |l0j which we, however, do not consider in the following. 
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B. Calculations 

1. Principle 

In the following, we use the structure equations ||HJ| and Q for our numerical calculations of the mass and the radius 
of white dwarfs and neutron stars. We need two initial values for the calculation, p(0) and m(0). From equation (J5J 
it is clear that m(0) = 0. We will increase r in small steps and for each step in radius the new pressure and new mass 
are calculated via: 

. . . Ge(r)m(r) . , . 
p(r + Ar) = c2? , 2 Ar+p(r) (55) 

47rr 2 e(r) 

m[r + Ar) = - Ar + m(r) (56) 

with Ar as a fixed constant which determines the accuracy of the calculation. This is the so called Newton method, 
whereas later we also employed a four-step Runge-Kutta method which allows a more sophisticated numerical inte- 
gration. As mentioned before in section I'll A II the pressure will decrease with increasing r. At a certain point, the 
pressure p(R) will become zero (or even negative) which will terminate the integration and determines the radius of 
the star R and as its mass M = m(R). For our computation we worked with Fortran and Matlab. 

2. Polytropic EoS 

First, we will use the polytropic equation of state 

p = Ke 1 . (57) 

In this subsection we closely follow p|. We will just briefly discuss the method which was used in as well as our 
results and refer the reader to the original paper for more detailed information. The complication in solving the 
differential eqs. JSJ and Q is to find a way to get for a certain pressure p(r) the corresponding energy density e(r). 
The polytropic form of the EoS, eq. (|57|l . avoids this problem and being plugged in eq. J3J gives: 

dpjr) _ -R p(r) 1/7 m(r) 
dr ~ r 2 K l / 1 ' [ ' 

with: 

m(r) GM Q 

m ( r ) = -in—, Ro = — — (59) 
M Q c 2 

where M Q is the mass of the sun, hence m(r) is dimensionless, and Rq is half of the Schwarzschild radius and is given 
in km. We also get a new expression for the differential equation for the dimensionless mass: 



dm(r) Airr 2 ( p(r) ^ 1 



dr Mqc 2 \ K 



(60) 



The units for energy density e(r) and the pressure p(r) are dyne/cm 2 , the mass is expressed in Mq and the radius in 
km. We can now calculate the mass and the radius of a white dwarf for a given central pressure. 

We will distinguish between the relativistic case and the non-relativistic case by setting 7 = 4/3 and 7 = 5/3, 
respectively. Choosing for the central pressure po = 2.33002 • 10 22 dyne/cm 2 , where this value satisfies the condition 
kp <C m e c (cf. eq. I |22| l). our numerical calculation with the Newtonian Code gives us the curve shown in Fig. ^ 
for the pressure-radius and the mass-radius relations. The white dwarf's mass is M = 0.40362 M Q and the radius 
R = 10,919 km. 

To satisfy the condition Uf S> m e c for a relativistic white dwarf, we have to choose a larger central density compared 
to the non-relativistic case, p = 5.619 • 10 25 dyne/cm 2 is an appropriate value (see Fig. 0). 

There is one interesting point as far as the relativistic polytrope is concerned. As we have seen in section III A 31 
the mass limit to a white dwarf must be 1.4312 M Q if we set rj equal 2 and use a relativistic polytropic form for 
our equation of state. As these are exactly the same underlying assumptions as in our numerical computation, the 
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i 1 1 1 1 1 1 1 r 




r in km 

FIG. 1: The mass m(r) and the pressure p(r) for a non-relativistic white dwarf with a central pressure of po = 2.2 ■ 10 
dyne/cm 2 . 




r in km 

FIG. 2: The mass m(r) and the pressure p(r) for a relativistic white dwarf with a central pressure of po = 5.62 • 10 25 dyne/cm 2 . 

maximum calculated mass must necessarily be the Chandrasekhar mass. Indeed, our numerical calculation gives us a 
white dwarf with a radius of 5710 km and a mass of 1.43145 M Q . With the chosen high initial pressure value we hit 
the Chandrasekhar mass limit very well. 

Note, however, that according to the expression l|51|) the mass does not depend on the central pressure. To see this 
in more detail we have employed different central pressures: Fig. 0] shows the calculated masses of white dwarfs for 
initial pressures from 1 • 10 23 dyne/cm 2 to 2.5 • 10 24 dyne/cm 2 . For this range of central pressures we always compute 
the same mass of 1.43145 M Q . Hence, the numerical calculations agree with theory. As an additional check we also 
plot the corresponding radii which we computed numerically (points referred to as "R") and calculated analytically 
with eq. I|48|) (referred to as "Rfit")- The deviations of the numerically computed radii from the analytical ones 
are about 0.2%. Taking a look at the plot, we find that for increasing central pressure the mass of the relativistic 
white dwarf stays constant while the corresponding radius decreases. The behaviour of the mass and radius of a 
non-relativistic white dwarf with smaller central pressures is shown in Fig. [21 for comparison. In contrast to the 
relativistic case the mass increases now with increasing central pressures. 

Silbar and Reddy already pointed out in [| that the two polytropic white dwarf models examined here are not quite 
physical, as the non-relativistic EoS just works for central pressures below a certain threshold while the relativistic 
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FIG. 3: Calculated masses and radii for white dwarfs using the non-relativistic polytrope. 
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FIG. 4: Calculated masses and radii for white dwarfs using the relativistic polytropic EoS. 



EoS is not applicable for small pressure. They suggest finding an EoS which covers the whole range of pressures by 
describing the energy density as a combination of a relativistic and a non-relativistic polytrope. We decided to try 
another way. On the other hand, we will employ not a guessed combination but a consistent relation p(e) . 



3. Looking for zero-points 

We remind the reader that it is our goal to find a relation between the energy density and the pressure which is 
needed for the calculation of the mass and the radius of the star. In the following, we want to use the full equations 
for a degenerate electron gas to derive the EoS consistently for the whole pressure range considered for white dwarfs. 
As we see from eqs. (|18fl and <|22|) . it is impossible to write explicitly the energy density as a function of the pressure 
without using certain approximations on e or p just like it was done in the previous section (i.e. the polytropic 
approximation) . 

On the other hand, a numerical solution to this problem is straight forward. As both the pressure and the energy 
density are functions of the Fermi momentum, we just have to find for a given pressure the corresponding Fermi 
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momentum and use it in the equation for the energy density. With this procedure we are able to calculate the energy 
density e for a given pressure p. To get the Fermi momentum we apply a root-finding method for the equation 

P (k F )-p = 0, (61) 

where p is given and p(kp) taken from eq. I|22|l. As long as we have a simple expression for the derivative of pressure 
we use the Newton-Raphson-Method. In other cases, when we have more complicated equation of states, we take a 
simple bisection method, which is less accurate, but easier to handle. In addition, we write the energy density as the 
sum of the nucleon energy density and the electron energy density (in the polytropic approximation the energy density 
of the electrons was ignored). Employing this numerical prescription to the solution for the EoS, we can calculate the 
mass/pressure and the radius/pressure dependence, just like in the polytrope case. 
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FIG. 5: Calculated masses and radii for white dwarfs using a consistent EoS. 

In section III A 41 we discussed the corrections from general relativity and we will now test their influence on the 
masses and radii of white dwarfs. Fig. [5] shows our calculations for three different cases. In the case referred to 
as "Polytrope" the energy density is given by just the nucleons, like in the polytropic EoS earlier. The "Newton" 
mass/pressure curve is calculated with an energy density consisting of both, the nucleon and the electron energy 
density (see eq. (116(1 ). for the "TOV" curve the EoS is like the one for "Newton" but with corrections from general 
relativity (eq. I|53[) ). All three curves confirm the Chandrasekhar mass limit. 

For pressures smaller than 10 27 dyne/cm 2 the three cases are indistinguishable while for larger pressures the curves 
start to deviate significantly from each other. We have indicated two points on the mass/pressure curves: the first 
one corresponds to a central pressure where kp — m e c, i.e. the point where the electrons start to become relativistic, 
the second point is at kp = 10m e c and marks a central pressure where the electrons are highly relativistic. At the 
latter point, one can already see how the three curves start to deviate from each other. The similar behaviour for low 
pressure is due to the fact, that the electron contribution to the energy density as well as the corrections from general 
relativity are negligible. This is the non-relativistic region and the behaviour of the curves is comparable to the one 
of the non-relativistic polytrope (p cx e 5 / 3 ), the mass increases with the central pressure. 

The difference of the three curves in Fig. [5] starts to appear around p — 10 27 dyne/cm 2 where kp — 10m e c. 
Nevertheless, in all three cases the masses stay almost constant at a value of about 1.4 Mq around po = 10 28 dyne/cm 2 . 
At this pressure the relativistic electrons dominate the mass/pressure relation like in the relativistic polytrope case 
(p cx e 4 / 3 ). 

We will discuss the behaviour at larger pressure starting with the case "Polytrope". As the energy density for this 
model stems solely from the nucleons, the results for mass and radius for higher pressures have to merge into that of 
the relativistic polytrope. This is confirmed by the constant mass of 1.4314 solar masses (equal to the Chandrasekhar 
mass) for central densities above po — 10 31 dyne/cm 2 . The mass in the center of the white dwarfs increases with 
larger central densities while the radius decreases, such that the total mass stays constant. 

For high pressures the contribution of the relativistic electrons to the energy density becomes important which is 
taken into account in the "Newton" case. The EoS for electrons only can be approximated by an ultra-relativistic 
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polytrope (see eq. (G2}): 

P = l -e . (62) 

With increasing central pressure the range in radius of the white dwarf where the electrons obey the ultra-relativistic 
limit enlarges. The deviation from the "polytrope" -curve is an effect of the growing fraction of ultra-relativistic 
electrons. The masses for high pressures are now smaller than the ones in the "Polytrope" case for the following 
reason: the degeneracy pressure of electrons for the same energy density will be lower in the ultra-relativistic polytrope 
compared to the "Polytrope" case due to additional contribution of electrons to the energy density. Hence, the star 
will support less mass and the mass decreases as a function of pressure beyond the maximum mass. For the "TOV" 
mass/pressure curve we take the former case adding effects from general relativity (see eq. I|53[) h At a pressure of 
Po ~ 10 28 dyne/cm 2 the mass curve of white dwarfs exhibits a clear maximum with a maximum mass being a little bit 
lower than for the Newtonian cases i.e. M max = 1.39 M©. The slightly reduced maximum mass is in accord with the 
fact that the corrections from general relativity strengthen the gravitational force. With increasing central pressures 
and smaller radii, corrections become larger. Therefore, the mass of the white dwarf decreases with increasing central 
pressures up to a pressure of about po = 2 • 10 37 dyne/cm 2 , here the mass starts to rise to a maximum value of about 
0.51 solar masses and then decreases again. The general relativistic corrections applied in the "TOV"-curve enhance 
the gravitational force, so that less mass is needed to create the given central pressure. 




FIG. 6: The mass-radius relation for white dwarfs with corrections for general relativity. 



Fig. shows the calculated masses for the "TOV" case in dependence of the radius. The central pressure increases 
along the curve starting from the right to the left in the plot. The curve begins at large radii and small masses and 
continues to smaller radii and larger masses. Passing the maximum mass the curve starts to decrease with radius. 
For higher central pressures the curve starts to curl. This behaviour is typical for large central pressures, seen also 
in many other mass-radius relations for different equations of state. A detailed analysis of this curled shape can be 
found in [HI and [T^. 

In the next subchapter we will discuss that white dwarfs of the "Newton" and "TOV" case which are beyond the 
maximum mass limit (i.e. for p 3> po) are unstable. It can be shown that polytropic equations of state with an 
exponent less than 4/3 cause an instability (see e.g. 0): an increase in energy density causes an insufficient increase 
in pressure to stabilise the star. 



4- Stability 

First we argue that white dwarfs are unstable if dM/dpo < following F° r our argumentation we employ the 
mass curve in Fig. [7| and reason at first why white dwarfs with po < p™ ax — po(M = M max ) are stable. 
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p in dyne/cm 2 



FIG. 7: Calculated masses and radii for white dwarfs with corrections for general relativity 



Let us assume that a white dwarf oscillates radially and therefore increases his central pressure to p' = po + Spo with 
dp > 0. In Fig. □ the mass curve shows white dwarfs in hydrostatic equilibrium i.e. the degeneracy pressure and the 
gravitational force compensate each other. If for po < p™ ax the central pressure is increased, the white dwarfs are only 
in hydrostatic equilibrium if their masses are larger than before. Since this is not the case, the surplus in degeneracy 
pressure over the gravitational attraction will expand the star back to his former state. For pressures p' — po — Sp a 
Fig. shows that the mass M(p' ) of a white dwarf in hydrostatic equilibrium is smaller than M(p ). That means 
that if the star expands due to radial oscillations the resulting lower central pressure can only be maintained if the 
gravitational force, i.e. the mass, decreases. As the mass of the white dwarf does not change the star has to shrink 
again to its former state. Therefore, white dwarfs with po < p™ ax and dM/dpo > are stable. 

Now we discuss the situation for a white dwarf with po > p™ ax and dM/dpo < 0. If the central pressure is lowered 
to p' = po — Spo then M(p' Q ) should be larger than M(p ) to bring the star in hydrostatic equilibrium (see Fig.0). So 
the gravitational force is not large enough to stop the expansion and the white dwarf is unstable to radial oscillations. 
For an increased central pressure p' = po + Spo the mass curve shows that stability is only given for a white dwarf with 
Po < p nax . As the mass of the white dwarf is larger than the one required for hydrostatic equilibrium, gravitation 
will cause the white dwarfs to shrink and finally to collapse. 

In [T^ |. Q as well as in Q, which we will follow now, a more formal approach is presented. Again, one considers 
radial oscillations of the white dwarf. The radial modes w„ can be derived by solving a Sturm-Liouville eigenvalue 
equation with the eigenvalues Lo n 2 . In doing so one gets a range of eigenvalues with 

uj 2 < loi 2 < oj 2 2 < ... (63) 

where oj n are the eigenfrequcncies to the radial modes and ujq corresponds to the fundamental radial mode. For 
the instability analysis the so called "Lagrangian displacement" £ (x,t) is introduced in ref. @ which connects fluid 
elements in the unperturbed state to corresponding elements in the perturbed state. The normal modes of an 
oscillation can then be written as 

f(x,<)=f(x)e^ , (64) 

where £ l are any Lagrangian displacements. Hence, the star will be stable for J 2 > and unstable for to 2 < 0. This 
means that for stability all modes need to have real eigenfrequencies. But because of l|63l) the stability of the star 
depends only on the sign of uiq 2 . 

The instability condition can be connected to the total Mass of any single star in in hydrostatic equilibrium (i.e. 
every star lying on our M(R) curve). The variation SM of the mass out of the state of equilibrium can be expressed 
as follows: 

M = M a + SM + 6 2 M , (65) 

where M is the total mass at equilibrium. For all stars in hydrostatic equilibrium it holds that the variation 
SM — 0. It can be shown (see 0) that stability is given for S 2 M > with an onset of instability at S 2 M — 0. For 
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8 2 M < the star is unstable. In the mass/pressure relation one finds minima and maxima, where dM/dpo — 0. For 
these cxtrcma it can be deduced that also S 2 AI = 0. This means these cxtrcma are the critical points for stability. 
The latter condition relates to uo 2 — and therefore signals a change in the stability of the corresponding normal 
mode. 

In the case of multiple critical points it is also possible that an unstable configuration can change back to stability at 
the following extremum. Therefore one has to take into account the sign of the change in the radius R (i.e. dR/dpo < 
or dR/dpo > 0) in addition. If an equilibrium configuration of a star is now transformed from the low-pressure side 
of the critical point to the high-pressure side, one can find the following relations between the change of the radius R 
with respect to the central pressure po and the radial oscillation modes: 

• dR/dpo < 0: the square of the oscillation frequency of an even numbered mode changes the sign (w n 2 — > —uj n 2 ) 

• dR/dpo > 0: the square of the oscillation frequency of an odd numbered mode changes the sign. 

For the "TOV" curve in Fig. [S] there are three critical points. The first critical point appears at the maximum 
around po = 10 28 dyne/cm 2 . As can be seen in Fig. [7] dR/dpo is negative in this critical point, therefore an even 
mode is changing stability. Because of and the assumption that for low densities the ocsillation frequencies of 
all modes are real (i.e. uj, 2 > 0), ujo 2 should change its sign here. Therefore, for central pressures between the first 
critical point and the next one, the oscillation in the fundamental mode is unstable. 

At po ~ 10 37 ' 5 dyne/cm 2 , the second critical point is reached. Here dR/dpo > and oj\ 2 becomes negative. At the 
second maximum at pq ~ 10 39 dyne/cm 2 again the oscillation in an even mode must change the sign and as uo\ 2 is 
negative already, U2 2 becomes negative now. This stability analysis can also be applied to the case of the Newtonian 
white dwarfs as well as to the neutron stars. As mentioned before, the sequence can become stable again; then at the 
second critical point the fundamental mode changes its sign back and a stable sequence of compact stars starts again. 
This so called third family of compact stars appears for a strong first order phase transition to e.g. quark matter 

0HG1 

In summary white dwarfs become unstable as soon as their mass is above the Chandrasekhar mass limit. This is 
exactly what is happening in the case of a core collapse supernova explosion. It is quite interesting to see that such 
an intuitive reasoning already indicates what is happening in real nature. 



III. PURE NEUTRON STARS 
A. Non— relativistic polytropic case 

In this chapter we start looking at pure non-relativistic neutron stars, which are described by an EoS of a Fermi 
gas of neutrons. In the following we drop the subscript F for labelling the Fermi momentum due to practical reasons. 
With the neutron number density 

= Arg (66) 

and the mass density 

p = n n ■ m n (67) 
we get in analogy to the Fermi gas for electrons (see eqs. l(T7|) - (H2l) : 

e(x) = ^[(2a; 3 + a;)(l+x 2 ) 1 / 2 -sinh- 1 (x)] (68) 

P 

with 



(«) . ±U*-*.K + W + (69) 



7 4 5 

x=— , eo = -rrj ■ 70 



In the non-relativistic case eqs. I|68|) and (|69|l simplify to: 



e(x) ~ pc 2 = |x 3 , p(x)^^x 5 . (71) 
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Hence, the EoS can be described by a polytrope: 

p(e) = K nrel e 5 / 3 (72) 

with 

h 2 ( Sir 2 \ 5/3 2fi cm 2 , . 

K nrel = — 2 = 6.428 • 10~ 26 — m . 73 

We start our numerical calculation with our Newtonian code and the results are plotted in Fig. |H1 We arrive at 
neutron star masses of around 0.5 Mq and radii of 20 km to 30 km. The general relativistic effects for neutron stars 
should be larger than for the white dwarfs. In Fig. [S] we compare the Newtonian with the TOV-code results. We see 
that the general relativistic corrections are at the order of a few percent and become more and more pronounced for 
more massive neutron stars. 




p in dyne/cm 2 



FIG. 8: The mass M and radius R of neutron stars in the non-relativistic case as a function of the central pressure po- The 
results from the TOV equation are compared to the Newtonian limit. 



B. Ultra— relativistic polytropic case 

For the ultra-relativistic limit, x 1, one finds for the energy density and for the pressure of nucleons: 

^) = f* 4 , PW = ^ 4 ^ p = e/3 (74) 

which is similar to the case of ultra-relativistic electrons discussed before. However, the numerical integration of the 
TOV equations fails here: The condition for terminating the integration (p < 0) can not be reached. The pressure 
converges monotonically as a function of the radius against zero, so that the radius of the star approaches infinity. 
The mass does not go to infinity because dm/dr also reaches zero asymptotically. The problem with the relativistic 
case is, that its solution has some principal inconsistencies. Running from the initial pressure po to zero through 
the whole star, the pressure has to pass the region where neutrons become non-relativistic, necessarily. The pure 
relativistic approximation is not valid for the whole density range of a compact star. 



C. Full relativistic case 



It is essentially to find an equation of state which is valid over the whole range of pressure and energy density. In 
other words: Eqs. l|68|l and l|69|l should be used for the calculation. Again we use a root-finding routine, as we have 
no explicit equation of state of the formp(e) or e(p) any more. The results for R(po) and M(po) are shown in Fig. [5] 
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FIG. 9: Mass M and radius R as a function of the central pressure po for pure neutron stars. 
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FIG. 10: The mass M for pure neutron stars as a function of radius R. The maximum mass is M = 0.712 Mq at a radius of 
R = 9.14 km. 

The mass-radius relation M(R) is plotted in Fig. 1101 where the maximum mass can be easily recognised. The 
heaviest stable pure neutron star has a central pressure of po = 3.5 • 10 35 dyne /cm 2 , a mass of M = 0.712 M Q and a 
radius of R = 9.14 km as first calculated by Oppenheimer and Volkoff in 1939 [131 (they arrived at a maximum mass 
of M = 0.71 M and a radius of R = 9.5 km for comparison). 

IV. NEUTRON STARS WITH PROTONS AND ELECTRONS 

Because the neutron is unstable, a neutron star will not consist of neutrons only, but in addition a finite fraction 
of protons and electrons will appear in dense matter. The electrons and protons are produced by /3-decay: 



■p + e +v e . 



(75) 
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There exist reverse reactions also, so that an electron and a proton scatter into a neutron and a neutrino: 

p + e~ — ► n + v e . (76) 

As all reactions are in chemical equilibrium, for a cold neutron star, the following relation for the chemical potentials 
holds: 

Mn = Mp + Me- (77) 

Neutrinos do not have a finite chemical potential (fj, u — 0), as the neutrinos escape from the cold star without 
interactions. Due to stability reasons (see e.g. ^3), ^ is assumed that the whole neutron star is electrically neutral, 
hence: 

n p = n e k p = k e (78) 

and the Fermi momenta of protons and electrons must be equal. With the definition of the chemical potential for an 
ideal Fermi gas with zero temperature, 

»i(h) = ^ = (k 2 c 2 +m"c 4 ) 1/2 , i = n,p,e, (79) 
atii 

and using eqs. l(77|l and (|75)l . one arrives at the following relation: 

(k 2 n c 2 + mlc 4 ) 1 ' 2 - (k 2 v c 2 + m 2 ^) 1 / 2 ~ (k 2 p c 2 + m \c 4 fl 2 = 0. (80) 
This equation can be solved for the Fermi momenta of protons k p as a function of the Fermi momenta of neutrons k n : 

t , _ \{k 2 n c 2 + m 2 c 4 - m 2 e c 4 ) 2 ~ 2m 2 p c 4 {klc 2 + m 2 c 4 + rnjc 4 ) + m 4 ^} 1 ' 2 

plnj_ 2c(fc 2 c 2 + rule*) 1 / 2 ■ [ ' 

The total energy and the pressure are now the sum of the energy and pressure of electrons, protons and neutrons, 
respectively: 

etot= 5^ £i, Ptot= 5^ Pi, (82) 
where 

ei ( ki ) = t(k 2 c 2 +m 2 c A Yl 2 k 2 dk (83) 

Pi (ki) = l^^ 3 J** ( k2c2 + m * c4 r 1/2k4dk - ( 84 ) 

For the case when there are no neutrons present, i.e. k n = 0, eq. I|81|l gives: 

kJO) = — — 1 - " — = 1.264 • 10- 3 m„c. (85) 

2m n c 6 

This means, that eq. (|81|l can only be used, if k p > 1.264 • 10~ 3 m„c, which is fulfilled for p > p cr j t = 3.038 • 
10 24 dyne/cm 2 . Below this pressure, no neutrons are present any more, so the "neutron star" material consists of 
protons and electrons only. Fig. 1111 shows the equation of state for a free gas of protons, electrons and neutrons 
in /3-equilibrium in the form e(p). Around the critical pressure, a notable kink appears in the equation of state, 
which can be identified with the strong onset of neutrons appearing in matter with their corresponding additional 



2— n,p,e 



contribution to the energy density. (One can show that 



dp 2 



= co, but ^ 



is continuous and thus the physical 



situation is not that of a true thermodynamic phase transition.) 

Using the new equation of state, we start our TOV calculation and compute masses and radii which are depicted in 
Fig. 1121 and Fig.^] respectively. The result seems to be surprising at a first glance but can be explained by looking 
at the proton to nucleon fraction integrated over every single star as displayed in Fig. 1131 

Below the critical pressure of p = 3.038 • 10 24 dyne/cm 2 , the neutron star consists of protons and electrons only by 
definition. Above this pressure, neutrons appear, resulting in a sharp bend of M(po) and R(po). Between the critical 
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FIG. 11: The energy density e as a function of the pressure p for neutron stars including protons and electrons. At p 
3 • 10 24 dyne/cm 2 neutrons begin to appear. 
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FIG. 12: Mass M and radius R as a function of the central pressure po for neutron stars including protons and electrons. 




pressure and p ~ 10 29 dyne/cm 2 , the stars have a tiny neutron fraction in the core, which increases with the pressure. 
The increasing neutron fraction is responsible for the decreasing masses. On the microscopic scale the relativistic 
electrons combine with the protons to form neutrons, thereby reducing the electron degeneracy pressure. For the 
same energy density the pressure will be lower, if neutron production is included, so only a smaller mass is supported 
by the star. At p ~ 10 30 dyne/cm 2 the proton fraction drops to n p /n to t < 10~ 3 and the star becomes an almost pure 
neutron star (hence the name). The proton fraction is so small, that it can be safely neglected. As a consequence the 
stellar objects become more compact, with a smaller radius and mass. For larger pressures (p > 10 30 dyne/cm 2 ), the 
mass increases again with rising pressure. The neutrons start to become degenerate, the star is stable because of the 
neutron degeneracy pressure. 

So the behaviour of M(p ) is understood through the change of the proton fraction n p /n to t(po)- But what about 
the unusual peak of R(po)7 Similar to the ultra-relativistic description of a pure neutron star, the situation at 
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FIG. 13: The proton to nucleon fraction n p /ntot and the neutron to nucleon fraction n„/n to t, integrated over every single star, 
versus the central pressure po. 
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FIG. 14: The Mass M as a function of the radius 7? for neutron stars including protons and electrons. The maximum mass 
star has a mass of M — 0.699 Mg, the corresponding radius is R = 9.24 km. 



po ~ 10 29 dyne/cm 2 is, that the pressure approaches zero very slowly, but never reaches it or falls below. Fig. 1151 
shows a star with po = 2 ■ 10 29 dyne/cm 2 , where this behaviour can be clearly seen. So it is necessary to make a 
cut-off (for this case at p = 10 19 dyne/cm 2 ), so that the integration terminates as soon as this pressure is reached. 
During the numerical computation the mass converges to M, while the radius grows bigger and bigger. This neutron 
star consists of a dense central core and a huge low-density crust. The protons and electrons are responsible for these 
huge radii: the energy density for the same pressure is much lower, if no neutrons are present (see Fig. Illfl . Because 
dp/dr (eq. (|53|l ) is a function of the energy density e, the pressure p(r) inside a star will decrease much slower, if there 
are only protons and electrons left. For more massive stars this crust disappears and the radius is determined by the 
neutron core. After all, such a huge crust does not seem to be very realistic for neutron stars. 
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FIG. 15: The pressure p along the radius r for a single neutron star including protons and electrons at po = 2 • 10 29 dyne/cm 2 



Note that for the most compact neutron stars (i.e. for po > 10 dyne/cm 2 ) protons appear again. The proton 
fraction approaches a certain constant limit at high densities. Using eq. (|81ll in the ultra-relativistic limit, k n c 3> m n c 2 , 
results in: 



kpikn) 



l/2k„ 



l/8n r , 



(86) 



n t ot 



(87) 



The proton fraction for large po in Fig. ^]is of course smaller than 1/9, as only the central core obeys the relativistic 
limit, and no protons are present for the regions in the star with smaller pressure. The maximum mass star has a 
proton fraction of n p /n to t — 0.02. 

Another important point is, that all stars in the region of 10 24 dyne/cm 2 < po < 10 29 dyne/cm 2 are unstable, 
which can be seen from the M versus po plot in Fig. El Hence, there are only two different kinds of stable objects: 
first the "hydrogen white dwarfs", like stars consisting of protons and electrons, up to M = 4 Mq and radii of about 
10 4 km, then the almost pure neutron stars with a maximum mass of M — 0.699 M and a radius R — 9.24 km at a 
central pressure po = 3.412 • 10 35 dyne/cm 2 . Compared to M — 0.712 M Q and R = 9.14 km for the pure neutron star 
of the last section, the differences are quite small. The minimum neutron star mass is found to be about 0.03Mq with 
a radius of around 1000 km; the radius rapidly changes for this mass range and quickly shrinks to values of 200-300 
km for slightly larger masses. We note that these values for the minimum mass of neutron stars are not far from the 
result of M m i n = 0.09M Q found in more elaborate calculations (see e.g. 0). 



MODELS FOR THE NUCLEAR INTERACTIONS 



A. Empirical interaction 



In this and the following section we use MeV and fm for the energy and distance units. Additionally c is set to 1 
and he = 197.3 MeV fm. 

Following we introduce an empirical interaction for symmetric matter (n„ = n p ) of the form (see also [Tsj'k 

= m N + (E Q ) u 2/3 + -u + -u a , u = n/no . (88) 

n 2 a + 1 

The dimensionless quantity a, and the quantities A and B with dimension of MeV, are the fit parameters of the 
model. In fact, only the third and the fourth term of eq. Q88|) describe interaction terms, (Eq) is the average kinetic 
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energy per nucleon of symmetric matter in its ground state: 



3 k 2 
(E ) = 6 L ° 



5 2mjv 5 2rriN 



3 1 {3ir 2 h 3 n 



2/3 



(89) 



The factor two in the denominator of the parenthesis comes in as an additional degeneration factor for the isospin 
(compare to eq. I|13jl), 

In the model the following values should be reproduced by the ansatz of eq. (|88|) : 

• the equilibrium number density no = 0.16 nucleons/fm 3 =>• fco = 263 MeV 

• the average binding energy per nucleon BE = (E/A — mjv)|„ =fl = (e(n)/n — mjv) \ n =n = —16 MeV 



_ n dp{n) 
dn 



400 MeV 



• the nuclear compressibility in nuclear matter K{tiq) = 9- 

After some algebra one gets the following parameters: 

A = -118.2 MeV, B — 65.39 MeV, a = 2.112, {Eq) = 22 A MeV. 
For known e(n) it is possible to determine p(n) via: 



2 it, \ 5/3 A 2 B(J 

- (E ) ^° /3 + ^ru 2 + — — i 

3 2 a + 1 



.(T+l 



(90) 



(91) 



We proceed now to extrapolate to asymmetric nuclear matter. 

First, the neutron to proton ratio is expressed through a parameter a: 



Tin 



N - Z 



a 



n A 

Now the kinetic energy is given by the sum of the kinetic energies of protons and neutrons: 



e K E{n,a) 



3 h 2 



3 k 2 



5 2m 



N 



5 2m,N 

,5/3 



= n (Ep) — (l + a) 0/J + (l-a) 



5/3 



(Ef) = ] 



3n 2 H 3 n^ 2/3 



(92) 

(93) 
(94) 

(95) 
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The difference in the kinetic energy density relative to the symmetric case can be expressed as a function of n and a: 



Ae K E (n,a) = e KE (n, a) - t KE (n, 0) 



(Ef){^ (l + af /3 + (l-a) 



5/3 



- 1 



This gives the following Taylor-Approximation: 



5 9 /_ a 2 



Ae KE (n,a) =n(E F ) ^a 2 ^1 , } _ 

This motivates that it is sufficient to keep corrections of order a 2 : 

E(n,a) E(n,0) 2 

— j — = — j — + a S{n). 

The function S(u), which describes the symmetry energy, is assumed to have the following form: 

S(u) = (2 2 / 3 - 1) (E ) (u 2 / 3 - F(u)) + S F(u). 



(96) 



(97) 



(98) 



(99) 



Sq is an input parameter (the asymmetry energy) , which describes the energy difference between pure neutron matter 
and normal symmetric nuclear matter at ground state density no- I n the following calculations, So is set to 30 MeV. 



22 



The only restrictions to the function F(u) are f(l) = 1, and F(0) = 0, in order to satisfy 5(1) = So, and 5(0) = 0, 
respectively. Besides this, F(u) can be chosen freely. We make the straight forward choice, 



F(u) = 



(100) 



The factor (2 2 / 3 - 1) (E ) can be explained as follows: For pure neutron matter, a = 1, the kinetic energy is not given 
by eq. any more, the factor 2 in the denominator of the parenthesis has to drop out. The factor (2 2 / 3 — 1) (E ) 
takes this into account: 

^1 = m N + (E ) u 2 ' 3 + K + JL-u° + (2 2 / 3 - 1) (E ) (> 3 - F(u)) + S F(u) 

= m N + (E ) u 2 ' 3 + (2 2 / 3 - 1) (E ) u 2 ' 3 + ±u + —u° + (s - (2 2 / 3 - 1) (E )) F(u) 

Z a + 1 v / 

= m N + 2 2 / 3 (E ) u 2 ' 3 + ^u+ —u° + (So - (2 2 / 3 - 1) (E )) F(u) (101) 

Z a + 1 \ / 

Thus the ansatz of eq. I|99|l implies the correct expression for the kinetic energy of pure neutron matter. Also the 
different contributions to the potential energy are separated from each other. 
With e(n, a) = nE(n,a)/A it is straightforward to calculate p(n,a): 



p(n,a) = n — 
an 



e(n,a) 



p(n,0) + n Q a 2 [(2 2 / 3 - 1) (E ) (2/3w 5/3 - u 2 ) + S Q u 2 



(102) 



As seen in the previous section, the proton fraction for stable compact stars is small. Hence, the parameter a is set 
to be equal 1 and a pure neutron star is considered in the following. 
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FIG. 16: The mass M as a function of the radius R for the empirical nucleon-micleon interaction. The maximum mass star 
has a mass of M = 2.792 Mq, the corresponding radius is R — 12.46 km. 



With p(n) and e(n) given, one can start again the numerical computation. The results are shown in Fig. 1161 The 
maximum mass star has a mass of M — 2.792 and a radius of R = 12.46 km. The neutron star is much more 
massive, if nucleon-nucleon interactions are included. As the equation of state is much harder, which means that the 
pressure increases faster with the energy density, the star can support a much larger mass. 

Another interesting fact is the behaviour of the mass for decreasing central pressure. In a certain region the mass- 
radius curve drops to smaller radii as the pressure and the mass decrease, until a radius of R ~ 9.5 km and a mass 
of M ~ 0.1 Mq is reached at a central pressure of p ~ 1.7 • 10 33 dyne/cm 2 . Passing this point, the radius starts to 
increase again. This special shape is due to the equation of state employed here. In Fig. El one sees that the EoS 
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FIG. 17: The equation of state for the empirical nucleon-nucleon interaction compared to the one of a free neutron gas. 



differs from the ones before: There is a smooth bump at p ~ 4 ■ 10 31 dyne/cm 2 . For pressures below this value an 
increase in the energy density leads to a very small increase in pressure, which is due to the attractive part of the 
interaction for matter close to the ground state density. In this pressure region a higher central pressure can only be 
generated by a smaller radius at a nearly unchanged mass. The slope of the EoS for pressures larger than the value of 
p ~ 4 ■ 10 31 dyne/cm 2 changes noticeable, so that also the mass-radius relation for the corresponding neutron stars is 
affected. For larger pressures the pressure increases much faster with the energy density, which results in the abrupt 
rise of the masses (see Fig. I16f) . The repulsive force of the interaction becomes dominant for nuclear densities above 
ground state density. Now an increase in energy density leads to a strong increase in pressure also and a larger mass 
is supported by the star. 

In fact, the EoS used here is not valid for small pressures and energy densities. Below nuclear matter density no it 
is not realistic to use a uniform matter distribution anymore. The formation of a lattice of nuclei and the occurance of 
certain possible mixed phases become important in the outer layers, the so called crust of neutron stars. For pressures 
below 10 33 dyne/cm 2 one has to use a special equation of state to describe the matter in the crust. Taking into 
account a realistic EoS for the crust would wash out the bump seen in Fig. El an d also the characteristic features in 
the mass-radius relation described above. A good introduction to neutrons star crusts can be found in e.g. 20] and 



B. Skyrme parameterised interaction 

In the following, the properties of nuclear matter are considered by applying the Hartree-Fock-method on a phe- 
nomenological nucleon-nucleon interaction taken from |22j (for further information see also |23|'). 
Consider a system consisting of N p protons and N n neutrons, which are described by plain waves: 

- -Le^xsX/- (103) 
V V 

Here, \i is t ne degree of freedom for the isospin, \s f° r the spin. The self-consistent Hartree-Fock equations read: 

e a <p a (x) =ttp a (x) + ^2 J d3 wpl(y) v ( x ,y) [ ( Pa'(y)fa(x) - <p a (y)v a >(x)} . (104) 

a' 

where ^2 a , is the sum over all occupied states. The potential is spin- and isospin-independent, so the sum can be 
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divided into a part for protons and one for neutrons: 

k ip) M 

E- = E- + E ( 105 ) 

a' k,s k.s 

As the sum splits up, one gets a separate Hartree-Fock-equation for neutrons and protons: 

k<?> 

e^^\x)=t^\x) + (2s + l)J2 I ' d*y[^)*{y)^{y)]V{x,y)^\x) 

k> J 

kg" 

- 1-E fdW k t(yW { k p) (y)}V( x ,y)^\ x ) 
k' J 

+ (2,s + l)]T / d 3 y[^*(y)^\y)]V(x,y)^\x). (106) 

by taking into account the spin-degrees of freedom g s = 2s + 1 . The neutron wave function results from replacing the 
superscript (p) by (n). 

Next, the interaction is introduced as a Skyrme-type parameterisation: 

V(x,y) = S 3 (x-y)(^t 3 n-t ), (107) 

where to describes the attractive two particle nuclear interaction, whereas t 3 n describes the repulsive (and density 
depended) many-body interaction which is the dominant one at high nuclear densities. Here, n denotes the total 
number density of the system: 

n=(N n + N p )/V = n n +n p . (108) 
The single-particle-energies for protons and neutrons can then be derived to be: 

4 P) - i§r + dtsn-t )±-(N p /2 + N n ) 



c 



(n) 



2m p v 6 'V 
k'i , 1 I 



^ k 2m n ' ^6"°'~ " K "V 

The total energy for a homogeneous system is given by: 



+ (-t 3 n-t )-(N n /2 + N p ). (109) 



N p N n 

^ = ^£(4 p) + *;) + ^£(4" ) + *0 (no) 



=► E BF = N p \3_ + N ^3_ + \ ( l tsn _ M _L(^ /2 + 27VpiV „ + Nl/2) _ (m) 

The next task is to determine the phenomenological parameters to an d t 3 . For n = n = 0.16 fm 3 the total energy 
shall have a minimum with Ehf/A = BE = —16 McV. From these conditions one arrives at: 

5m n c 2 \ 2 J 

16E B 16 (he) 2 /3tt 2 \ 2/3 . . m 
3n 15 TO„tr \ 2 / 

=>■ t = 1024.1 McV fm 3 , t 3 = 14600.8 MeV fm 6 . (113) 
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p in dyne/cm' 



FIG. 18: The equation of state for the Skyrme-parameterisation and for a free neutron gas. 



With 

dp{n) 



K(n ) = 9- 



dn 



(114) 



one finds a value of 376.4 MeV for the nuclear compressibility, which is less than for the empirical interaction. 
Again, we assume pure neutron matter, i.e. we set N p = and N n = A: 

Ehf 3 h„ 1 , 1 x . » 

-4 £ = i^f- + 7 «*3n-t n. 115 
A 5 2m n 4 6 

With the total energy density e(n), including the rest mass term m n • n, 

e(n) - (m n + E HF /A) -n (116) 

= mnn + _J_( 37r 2 ?i 3 ) 2 / 3 n 5 / 3 + ^n 3 ~ , (117) 
10m„ v y 24 4 ' v 7 

one can derive p(n) according to the thermodynamic relation: 

p(n) = u— e(n) - e(n) (118) 
an 

=> p(n) = -^(3.^ 3 ) 2 / 3 n 5 / 3 + ^n 3 - |n 2 . (119) 
10m„ 12 4 

The resulting equation of state is plotted in Fig. 1181 which appears to be very similar to the empirical interaction 
discussed before. 

Fig. E| depicts the corresponding mass-radius plot. The maximum mass star has a mass of M = 2.438 M 
and a radius of R — 10.98 km. These values are only a bit smaller than the ones calculated before for the empirical 
interaction, also the shape of the curve looks similar. The sharp turnover in the lower left corner of the plot corresponds 
to a minimal radius of only R ~ 8.6 km and a mass of M ~ 0.1 M Q . 



C. Comparing the two models 



It is instructive to take a closer look on the two different EoSs considered here, as the results do not differ that 
much from each other. The two energy densities (eqs. I|98|l and l|117|) l are: 



26 



3,0-. 
2,5- 
2,0- 

o 

.£ 1,5- 
1,0- 

0,5 
0,0 




10 



— r - 
11 



R in km 



12 



'T- 
IS 



— 1 

14 



FIG. 19: The mass M as a function of 7? for the Skyrme-parameterisation. The maximum mass star has a mass of M 
2.438 M and a radius of R = 10.98 km. 



• empirical nucleon-nucleon parameterisation: 

= mn + {Eo) u 2/3 + A + _B_ a + (2 2 /3 _ 1} {Eq) ( 2/3 _ \ + SqU 

n 2 a + 1 V / 

• Hartree-Fock method and Skyrme parameterisation: 

£HF{n) 3 2*-3n2/3 2/3 , 1 A , . \ 

= mn + iW 37r n ] n + i ( 6 t3?l ^ io)n 

With n = u • «o and eq. (|89|l one can rewrite the two equations in the following form: 

= mn + 2 2/3 ( £ Q> w 2/3 (2 2/3 _ 1} {Eo) + S \ u + B 

n \ 2 / cr + 1 

= TO „ + 2 2 /3 (£ ; o)w 2/3_ 1 nQU+ l ian 2 u 2_ 

71 4 24 

By inserting the numbers (|113|l . the energy density exhibits the following parametric form: 



>(n) 



n 

e HF (n) 



n 



= (939.6 + 35.1u 2/3 - 42.1m + 21.0u 2 ' 112 ) MeV ^ 21.0u 2 ' 112 MeV 
= (939.6 + 35.1u 2/3 - 41.0m + 15.6u 2 ) MeV ^ 15. 6u 2 MeV. 



(120) 
(121) 

(122) 
(123) 

(124) 
(125) 



If the energy density follows a power-law as a function of n, it is easy to find an explicit equation of state of the form 
p(e): 



c ■ n 



de a -,\ k 

p = n ■ — e=(fc — l)c-n 

an 

p(e) = (fc-l)e . 



(126) 



In our case, k = 3.112 for the empirical interaction and k = 3 for the Skyrme-parameterisation. One can see, that 
the first EoS is stiffer than the second. This explains why the masses in the empirical model are bigger as the masses 
in the Skyrme-parameterisation, the harder EoS can support more mass. 
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p in dyne/cm 2 

FIG. 20: The speed of sound depending on the central pressure for the two different models. The maximum mass stars are 
marked with a dot. 



D. Speed of sound and causality 



The speed of sound in our neutron star must be less than the speed of light to fulfil relativistic causality. The 

dp 

777 



adiabatic speed of sound is given by c s = ^ in a non-dispersive medium. Hence, we have: 



Using eq. I|12rj|) one gets: 



dp 



1 < 1. (128) 



Hence, the power coefficient k must be less than or equal two. So both EoSs will become acausal for large n. This 
behaviour is demonstrated in Fig. [201 where the speed of sound is plotted against the central pressure. As one can 
see, for the maximum mass star the speed of sound c s /c is larger than one for both interactions considered, hence, 
correction factors have to be added to the interaction, which include (or simulate) relativistic effects (see e.g. 0). 



VI. SUMMARY AND BRIEF OVERVIEW ON PRESENT DAY DEVELOPMENTS 

In chapter m we discussed the common structure equations for stars as well as the relativistic and non-relativistic 
equations of state (EoS) for a Fermi gas of electrons and nucleons. We simplified the EoS to a polytropic form and 
derived semi-analytically the Chandrasekhar mass limit of M = 1.4312 (2/rj) M for white dwarfs. We calculated 
numerically mass-radius relations for white dwarfs using the relativistic and non-relativistic polytropic approxima- 
tions. For the relativistic case we find the computed masses to be constant and in very good agreement with the 
Chandrasekhar mass limit. The numerical results were also checked by the analytical results for white dwarf radii. 

We introduced a numerical method to calculate the EoS of a free gas of particles for the general case, i.e. without 
using the polytropic approximation, and computed the corresponding masses and radii for white dwarfs. In this case, 
the calculated masses became larger with increasing central pressures up to a maximum mass of around 1.4 Mq. For 
high pressures, the masses stayed constant at the maximum value showing that the EoS turned into a relativistic 
polytrope. 

We calculated the mass-radius relations without as well as with corrections from general relativity. The mass curve 
showed an increasing mass for larger central pressures up to po ~ 10 28 dyne/cm 2 reaching a maximum mass of about 
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1.40 M© without and one of 1.39 M Q with relativistic corrections. The slightly reduced maximum mass in the latter 
case was in accord with the fact that general relativity corrections strengthen the gravitational force. 

We discussed also that white dwarfs and neutron stars are stable against radial oscillations for masses smaller than 
the maximum mass but unstable beyond the mass peak. 

In chapter ITTT1 we looked at pure neutron stars using an EoS of a free Fermi gas of neutrons. We computed the exact 
equation of state by a root finding routine. The maximum mass for such a neutron star was found to be M — 0.712 
M Q with a radius of R = 9.14 km. Then, we included protons and electrons in /3-equilibrium to the composition of 
neutron star matter. The maximum mass of the corresponding neutron star changed slightly to M = 0.699 M Q and 
the radius to R = 9.24 km. 

Next, we focused on the impact of nuclear interactions on the properties of pure neutron stars. We analysed 
two different types of interaction, an empirical interaction and the Hartree-Fock method applied to a Skyrme- 
parameterisation. As the two resulting EoSs were very close to each other, the results for the mass and radius curves 
were very similar, but quite different to the cases of free Fermi gases considered before. The maximum masses of 
neutron stars reach now values of M — 2AMq and M = 2.8Mq, respectively, which are up to four times larger than 
for the free cases. The corresponding radii increase to R = 11 to 12 km, respectively. The repulsion between nucleons 
stiffens the EoS, increases the pressure at fixed energy density, and enables to stabilise more massive neutron stars 
compared to the free Fermi gas case. Caveats of the simple nuclear interaction models used were exhibited, as the 
neglect of relativistic effects and the EoSs becoming acausal for large central pressures. 

We close this section with a brief discussion of present day developments in the field of compact stars: astrophysical 
data, the EoS of dense matter and its impact on the global properties of compact stars (for a recent general introduction 
on latest developments in the field we refer to H^). 

More than 1500 pulsars, rotating neutron stars with pulsed radio emission, are known today. About seven thermally 
emitting isolated neutron stars have been observed by the Hubble Space Telescope (HST), by the x-ray satellites 
Chandra and XMM-Newton and by ground based optical telescopes as the European Southern Observatory ESO. For 
the closed and best studied one, RX J1856. 5-3754, the measured spectra can be well described by a Planck curve (2||. 
As the outermost layer of a neutron star is the atmosphere up to a density of 10 4 gm/cm 3 consists of atoms, their 
presence should show up in the spectra. However, there is a lack of any atomic spectral line in the spectra which is 
not understood at present (27). For three selected rotating neutron stars (pulsars), nicknamed the Three Musketeers, 
one is even able to make phase resolved studies of the spectra, which allows to determine that there are hot spots on 
the surface of the neutron star [28( ■ 

The high-density EoS for neutron stars and the appearance of exotic matter has been studied in modern relativistic 
field-theoretical models fl3| . In particular, hyperons will appear around twice normal nuclear matter density (see 
e.g. [2£j and references therein). The possibility of Bose condensation of kaons in dense, cold matter has been 
investigated also |30| . Last but not least, the phase transition to strange quark matter, quark matter with strange 
quarks, will appear in the high-density limit |3lJ,|32j. If the phase transition of these exotic forms of matter is strongly 
first order, a new class of compact stars emerges as a new stable solution of the TOV equations. While the existence of 
this third family of compact stars besides white dwarfs and neutron stars was known earlier ll4 |. |3 3L |3 4l|. th e solution 
was rediscovered in modern models for the EoS and a Gibbs phase transition only recently |l5l Il6t l35l l36j . The 
topic of the third family of compact stars was a subject for another undergraduate student project of Jean Macher 
and one of the authors |37l |. The study of dense quark matter, as possibly present in the core of neutron stars, has 
been the object of intense investigation during the last few years due to the rediscovery of the phenomenon of colour 
superconductivity. Here, two colour charged quarks form a (colour) superconducting state and change the properties 
of quark matter (for reviews see e.g. [38l 133 . lioj ) . Neutron stars with colour su per conducting quark matter have been 
calculated by many authors in the last couple of years, see HJ, |H, |H, H3, HHH, H3, H El H3, IsjJ. The physics of 
strange quark matter in compact stars and its possible signals are reviewed most recently also in |3ll 132) . 

The physics of neutron star has also impacts not only on pulsars but also in other active research fields of astro- 
physics. For example, in x-ray binaries a neutron star accretes matter from his companion, which is an ordinary star 
or a white dwarf. When the matter falls onto the neutron star, x-ray bursts occur. Redshifted spectral lines have been 
measured from such x-ray binaries, which allow to constrain the compactness and by that the EoS of neutron stars 
|52| . Core collapse supernovae form a hot proto-neutron star first which cools subsequently and emits a neutrino wind. 
That neutrino wind influences the production of heavy elements in the supernova material in the so called r-process 
nucleosynthesis (see e.g. |53| and references therein). Neutron stars can collide with a corresponding emission of grav- 
itational waves, which will be measured by gravitational wave detectors in the near future, and gamma-rays, being 
considered a prime candidate for gamma-ray bursts and another site for r-process nucleosynthesis. The dynamics of 
the collision depends on the underlying EoS of neutron stars ..54] ! In the future, the square kilometre array (SKA) 
will measure more than 10.000 pulsars, with expected 100 binary systems with a neutron stars |55|. so the prospects 
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are bright for learning more about the physics of neutron stars! 
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